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Using the analytic and experimental techniques we present an exploratory study of the mass 
distribution features of the high coincidence of centre of mass of heterogeneous bone tissue in vivo 
and its centroid of geometry position. A geometric concept of the average distribution radius of 
bone issue is proposed and functional relation of this geometric distribution feature between the 
partition density and its relative tissue average distribution radius is observed. Based upon the 
mass distribution feature, our results suggest a relative distance assessment index between the 
center of mass of cortical bone and the bone center of mass and establish a bone strength equation. 
Analysing the data of human foot in vivo, we notice that the mass and geometric distribution 
laws have expanded the connotation of Wolff's law, which implies a leap towards the quantitative 
description of bone strength. We finally conclude that this will not only make a positive contribution 
to help assess osteoporosis, but will also provide guidance to exercise prescription to the osteoporosis 
patients. 
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Bone tissue structure and function are largely asso- 
ciated with its mechanical and biological environment 
[l|, [l, 0, 0, m. Growth, modeling and remodeling are 
the basic physiological features of bone. Biomechanically, 
bone growth is defined as mass changes @ , and mechan- 
ical force and movement play a role in bone growth 0. 
When the skeleton bears loads externally, bone tissue will 
undergo adaptive changes such as reabsorption or remod- 
eling [8| , and point-to-point changes to the material prop- 
erty by changing mass distribution 

SESIiil will maxi- 
mize its external loads. Individual bone growth indicates 
that external force has great effect on cross-sectional ge- 
ometry and internal anatomy [l2[. Bone structure is 
an optimization of stress transformation [l| and it is an 
adaptive response to incorporation of minimal weight to 
maximal strength by some special rules [13] • Bone phys- 
iological activity is regarded as a process of optimization 
0, [3, [H, Q • Consequently, stress has caused adaptive 
changes of bone shape and structure, which involve con- 
stant optimization of structures. But it remains unclear 
what distribution principle these changes follow. 

From the biomechanical perspective, osteoporosis 
means a sharp drop of bone mass and strength and they 
cannot meet the demands of adaptive strength and move- 
ment load [13] • Many mechanical models adequately rep- 
resent the relation between the bone geometry and its 
strength [H, 0, Q , as well as the correlation between 
bone density and its strength [2l| . While the phenomeno- 
logical models may often be helpful in obtaining a qual- 



itative understanding of the data, microscopically these 
models do not provide a trustworthy guide into unknown 
territory of an accurate relation between the distribution 
of bone tissue and its strength. Our approach will be to 
use a combination analysis of analytic and experimental 
techniques to examine the the above two uncertain ar- 
eas by setting up a bone strength equation to obtain a 
quantitative description macroscopically. 

Centroid of geometry (hereinafter referred to as COG) 
and the center of mass (hereinafter referred to as COM) 
of the homogeneous materials are coincide, whereas in 
most cases those of heterogeneous do not. Using CT 
(computed tomography) scan technology, we conducted 
the analysis to explore the relation between COM and 
COG of bone in the physiological activities, such as con- 
tinuous modeling and remodeling in its adaptive mechan- 
ical condition. At small enough CT resolution rate and 
its slice distance, the bone tissue density of infinitesimal 
bone volume segmentation, dV, can be regraded as con- 
tinuous. Its point density pi can be taken as that of the 
micro-element, drrii = pidV ^ thereby approximating the 
bone as a collection of particles. To find the optimal pro- 
gram for heterogeneous bone tissue mass distribution, we 
minimize 



min ^(pc) 



PiAV {xi-Xc) -^iVi-yc) -^{zi-ZcY 
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(1) 

where Pcixc^Vc^ Zc) and piixi^yi^ Zi) refer to the rela- 
tive locations of coordinators of bone COM and ran- 
dom point related to CT image, respectively. The se- 
ries Pi^y^ { \ - \ ^ \ Vi - Vc \ ^ \ Zi - Zc \ 

),^PiAl/( \ Xi - Xc \ ^ \ Vi - Vc \ ^ \ Zi - Zc \ ) ^re 
all convergent. Thus in the limit AV (AV = a6c. 
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TABLE I: Basic Information of the subjects. 



Sample size 

Age(year) 

Height (cm) 
Body mass{kg) 
Calcaneus volume(cm"^) 
Calcaneus density (g^/mZ) 



Wrestlers Volleyballers Seniors 
8 8 2 

21.00 ± 2.78 21.88 ± 0.99 64.50 ± 4.95 

168.00 ± 5.68 183.94 ± 3.90 150.50 ± 3.54 

65.52 ± 5.16 71.80 ± 5.20 52.88 ±3.15 

71.79 ± 7.86 81.79 ±4.26 49.43 ± 5.22 

1.47 ±0.04 1.49 ± 0.05 1.28 ±0.03 
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Volleyballer Senior females 



a ^ 0, 6 ^ and c ^ 0), the Abehan theorem on series 
Q leads to 

P = Pc. (2) 

where we have used p{x^y^z) for the COG. This shows 
that the optimal program for the heterogeneous bone tis- 
sue mass distribution should be the coincidence of its 
COM and COG. This signature does not, however, in- 
dicate a similar behaviour of the tissue distribution of 
bone in vive. In order to study the relation between the 
COM and COG of bone in vive, a CT scanning^ is con- 
ducted to the foot of eight volleyballers, eight classical 
wrestlers and two senior females. The position of COG 
is calculated by 

and that of COM by 

E^p E^p E^p 
Ep ' Ep ' Ep 

The comparative testing accurate distance between COM 
and COG positions is evaluated using 

^/{x- Xc)'^ + - Vc)'^ + - 

where a = b = y XY the CT image resolution, c 
the slice distance. 

Fig. [1] collects and displays our data shown in Table H 
The behaviour seen in Fig. [1] confirms our signature that 
the high coincidence of the position of COM and COG of 
heterogeneous bone in vive is independent of the chang- 
ing process of bone in its adaptive mechanical environ- 
ment and that the bone tissue mass distribution observes 
the optimal principle with a coincidence between COM 
and COG. 



^ The scanning is done by Philips/Brilliance 64 {120kv, Pixel Size: 
0.328 - 0.475mm, Slice Distance 0.330 - 0.450mm). 



FIG. 1: The position relation between COM and COG. The 
green ball refers to the COG of calcaneus while the red, black 
and orange balls represent positions of calcaneus COM from 
the wrestlers, volleyballers and the senior females (the colored 
ball radii are all ^\^a'^ -\- b'^ -\- c'^). 

The above signature brings us to the issue of mass 
distribution index of bone tissue. The mechanical prop- 
erties reveal that the elastic property and pressure- 
bearing strength of cortical bone is several times more 
than those of the same- volume spongial bone Us- 
ing the mechanical insight, we divide the tissue con- 
tinuous density into three parts; bone marrow, spon- 
gial bone and compact bone with their COM and COG 
highly coincident as indicated by the observed fact p = 
Pc- Using {xci - XcY + (^ci - Vc?' + (^ci - ^c)^, where 
[p^ci-, Vci-, Zci) refer to COM of bone tissue and (xc, yc^ ^c) 
refers to COM of calcaneus, we calculate the distance 
between each individual tissue COM and that of the cal- 
caneus. In order to avoid the effect from the size of the 
subject's calcaneus, we standardize the average distribu- 
tion radius of calcaneus so as to compare calcaneus of 
different volume, thereby establishing an anastz for the 
distribution index: 

E(\/(^j - + [vj - vcY + {zj - zcY) ' 

where i = 1,2,3 stands for bone marrow, spongial bone 
and compact bone respectively, j refers to all the tissues 
that make up bone and c stands for COM. 

From Fig. [2] it is clear that the position of compact 
bone COM of the volleyballers is closest to that of the 
calcaneus COM. Distinct difference exists between the 
volleyballers and wrestlers, so is the difference between 
the senior females and the volleyballers and wrestlers. 
The distance of the senior females is the largest. If such 
a trend continues for larger sample sizes, it will add one 
more quantitative evaluation index while diagnosing os- 
teoporosis. 

To see whether the changes of BMC (bone mineral 
content) and BMD (bone mineral density) bring about 
changes in the geometric distribution of bone tissue si- 
multaneously, we segment the density pi and use = 
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□ Wrestler □Volleyballer ^ «;*.n.of ♦♦•r^^* 

FIG. 2: Relative positions of bone tissue COM and the cal- 
caneus COM. 



o Wrestler □ Volleyballer a Senior females 
= 0.993 = 0.990 - 0.991 
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Bone tissue density 



FIG. 3: Relation between bone tissue density and its average 

distribution radius. Density is converted by p = PcrayVaiues 

^ PH20 

and distribution radius standardized by pi. 



to calculate the average distribution radius of the 
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relative tissue of segmented density to the calcaneus 
COM. When the continuity of segment density is guar- 
anteed and the calcaneous density variation ranges are 
fixed, the relation between the calcaneous density and 
the average bone tissue distribution radius has been de- 
veloped. After an analysis of the fitting function, we 
establish a functional relationship between the two for a 
typical value of correlation coefficient R > 0.99. Fig. [3] 
shows the obvious difference of average distribution radii 
between the senior females and the athletes. For the den- 
sity > 1.7, the average distribution radius of the senior fe- 
males drops dramatically. When comparing the wrestlers 
and the volleyballer s, when the density is greater than 
2.0, the bone tissue average distribution radius of the 
volleyballer s is more than that of the wrestlers. That is 
to say, when the geometrical distribution of bone tissue 
follows Eq. ([2j), senior females' bone loss will be accom- 
panied by a decrease of distribution radius of compact 
bone. When comparing the volleyballer s and wrestlers, 
there will be an apparent increase of distribution radius 
of compact bone. We may safely say that the bone adapt- 



ability includes not only changes of BMC and BMD, but 
also changes in radius of bone tissue distribution. Once 
the bone shape, tissue density and their corresponding 
volume have been determined, the tissue geometric dis- 
tribution will determine the strength of bone tissue. The 
moment of inertia of bone is an important index to re- 
flect bone strength. For the heterogeneous materials the 
bone density and intensity follow a non-linear relation- 

k 

ship. We introduce a coefficient and combine the seg- 
mented density and intensity to calculate the segmented 
strength on the basis of moment of inertia 

and establish a functional relation between the calcaneus 
bone strength and the segmented density as 



= f f{p)dp, 

J a 



(4) 



where a refers to strength of calcaneus, f{p) = 
pVr'^{p)e^ . When Ap is small enough, and f{a) ^ /(^), 
Eq. [4] can be approximated by 



Pmin) 
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pmax and pmin ^cfcr to the maximal value of compact 
bone density and the minimal value of spongial bone den- 
sity, respectively and exp(p^) is the coefficient parameter 
of bone density and strength. Note that the above equa- 
tion holds for continuous heterogeneous material only. 

Fig. m shows that when the bone tissue density of 
wrestlers and volleyballer s is greater than 1.8, the differ- 
ence in intensity between the two grows and reaches its 
maximum in the range of 2.4 — 2.5. On the other hand, 
the bone intensity of the senior females begins to show 
a larger difference with that of the athletes for density 
> 1.4. 

Concerning bone tissue distribution, our study con- 
firms that the high coincidence of COM and COG of 
heterogeneous material. This coincidence is the prereq- 
uisite to meet the requirement of max-min-principle and 
forms the bases for developing an evaluation index. We 
noticed that there is no difference in spongial bone be- 
tween the wrestlers and volleyballer s, whereas there is 
obvious difference in compact bone. This would explain 
the movement of the compact bone COM towards the 
calcaneus COM, which enables the calcaneus structure 
to bear greater stress. This can also be a representa- 
tion that bone can yield adaptive changes functionally. 
What's more significant is the fact that the compact bone 
COM of the senior females moves away from the calca- 
neus COM. If a larger size sample can verify this, it will 
bring greater significance to the clinical practice. 

One of the main aims of bone study is to conduct qual- 
itative analysis to bone strength. Bone strength relies on 
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o Wrestler = 0.994 

□ Volleyballer r' = 0.970 

A Senior females = 0.980 
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Bone tissue density 



FIG. 4: Geometric distribution of bone strength. The tis- 
sue volume of the segmented density has been normahzed 



by Vi 
fi = 
1.68. 



y , and its distribution radius standardized by 
The typical value of parameter k is chosen to be 



bone mass, but the uni-index of bone mass cannot paint a 
holistic and realistic picture of bone intensity objectively 
2qI , 21 , 2^. Modeling analysis reveals that the main fac- 
tors that determine the structure strength include mass 
distribution, geometric distribution and moment of in- 



ertia of various tissues [2J, |25|. Eq. ^ has success- 



fully combined those factors and from a mathematical 
perspective it has illustrated that volume, mass, density, 
distribution radius and moment of inertia of bone tissue 
cannot be employed individually to assess bone strength. 
Eq. (|4]) also reveals that criteria of selecting a training 
approach to increase or improve bone mass of compact 
bone and its distribution radius that are clinically signif- 
icant. The bone tissue geometric distribution principle 
sheds light on the effect to bone structure from differ- 
ent types of training. A geometrically-distributed bone 
strength equation can mirror the effects to bone strength 
from various factors. When the physiological bone mass 
decrease has become unavoidable, will exercise patterns 
be able to change its geometric distribution? If yes, Eq. 
(|4|) will undoubtedly provide some guidance to the de- 
velopment of exercise prescription and it can also be em- 
ployed as an important evidence to examine and modify 
exercise prescription. 
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